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Abstract: 

Gravitational contraction always generates a radially directed momentum flux. A 
particularly simple example occurs in the electron-degenerate cores of AGB stars, which 
contract steadily under the addition of helium ashes from shell hydrogen burning. The 
resulting momentum flux is quantified here. And since the cores of AGB stars lack efficient 
momentum cancellation mechanisms, they can maintain equilibrium by exporting their 
excess momentum flux to the stellar envelope, which disposes of much of it in a low 
velocity wind. Gravitational contraction easily accounts for the momentum flux in the solar 
wind, as well as the flux required to lift mass into the dust formation zone of every AGB 
star, whereon radiation pressure continues its ejection as a low velocity wind. This 
mechanism explains the dependence of the AGB mass-loss rate on core mass; its 
generalization to objects with angular momentum and/or strong magnetic fields suggests 
a novel explanation of why most planetary nebulae and proto planetary nebulae exhibit 
axial symmetry. 

Quasistatic contraction is inherently biased to the generation of the maximum possible 
momentum flux. Its formalism is therefore readily adapted to providing an upper limit to 
the momentum flux needed to sustain mass loss when this begins from a semi-continuous 
rather than impulsive process. 

Gravitation 
Stars: Mass-Loss 
Stars: AGB & post-AGB 
Sun: Solar Wind 
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I. Introduction 

Mass loss from the near environs of self gravitating bodies is ubiquitous. It occurs in 
spectacular fashion in the bipolar jets emanating from AGN such as M82, NGC4258 & 
M87, from quasars such as 3C273, as well as from stellar black holes such as SS433. 
Bipolar flows frequently occur in regions of star formation, as well as in the terminal stages 
of stellar evolution where they are found in most planetary and proto planetary nebulae. 
However, mundane, almost spherically symmetric winds emanate from many young star 
forming regions, as well as from massive main sequence stars, from the Sun, and from the 
AGB stars, which are commonly seen as large amplitude Mira variables. Yet our current 
understanding of the causality behind all of these mass-flows is slight. 

Mira variables lose mass to interstellar space in low velocity winds. In most stars these 
are driven by radiation pressure acting on dust grains which condense within their flows 
at a few stellar radii (Gilman 1972; Kwok 1975; Jura 1984). Direct observational evidence 
implicating dust in this process includes (i) a correlation of mass-loss rate, m, with UV 
opacity in carbon-rich circumstellar shells (CSs) (Knapp 1986); (ii) a strong increase in 
the expansion velocities, v^^p, of the winds from low mass Miras as their CSs become 
opaque, which facilitates the deposition of photon momentum (Lewis et al 1990); (iii) 
differences in the expansion velocities of carbon and oxygen rich stars that seem to 
depend only on envelope flavor (Lewis 1991); (iv) a difference in expansion velocities 
between OH/IR stars in the Galaxy and the LMC that are a result of their differing 
metallicities (Wood ef a/ 1992). The dense winds from OH/IR stars of luminosity L even 
carry a momentum flux, m v^^^ = x{L I c) with t»1: this has been shown by Netzer & 
Elitzur (1993), by Ivezevik & Elitzur (1995), and independently by Habing, Tignon & 
Helens (1994) to be explicable as a result of the multiple scattering of the stellar photon 
flux by dust particles. But 50 to 95% or so of the momentum flux, m v^^^, needed to lift 
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material from a Mira photosphere to infinity is consumed over the radial range below that 
at which dust forms, where the stellar atmosphere is essentially transparent. Indeed some 
AGB stars, such as % Cyg, have almost no dust, though they have winds. Consequently 
most of the momentum flux needed to support AGB mass-loss must come from within the 
star itself, and so has no direct dependence on the existence of dust. 

The detailed mechanisms for coupling momentum into a red giant's mass-loss stages 
remain uncertain (Hearn 1990), despite progress with the numerical simulation of mass- 
loss driven by Shockwaves in distended stellar atmospheres (Bowen 1988). Bowen's 
approach succeeds in pushing mass from the photosphere out to the cooler radii at which 
dust condenses: radiation pressure on the dust then sweeps the material away to infinity. 
But Bowen's model assumes the existence of a periodic piston below the base of the 
atmosphere to inject the requisite momentum flux: the existence of a flux-source within the 
star is a given of his models. Convincing ab initio calculations for the complete structure 
of AGB stars are still to come. In consequence many simple questions about Miras 
remain unanswered, such as what triggers their characteristic pulsation (Wood 1990). 
More fundamentally there is still confusion about exactly why a red giant star is a giant 
(Applegate 1988, Whitworth 1989). However, none of the past work addresses the 
question of the physical origin of the momentum-flux carried in the wind from the photo- 
sphere, which can be likened to (and may be identical to) the absence of an accepted 
physical causality for the momentum flux carried by the solar wind (e.g. Holzer 1976). 

In seeking an explanation for this flux I am drawn to the correlation between the 
occurrence of active phases of gravitational contraction in stars and the incidence of 
significant winds. Spherically-symmetric, 5-30 km s"^ winds occur during an early stage of 
protostellar evolution as well as the later stages of red giant evolution. This focus on 
contraction is relevant to AGB stars, since they occupy the stellar evolution niche arising 
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when steady contraction of the electron-degenerate core occurs, as helium ashes from 
shell-hydrogen burning make it ever more massive. Degenerate cores must contract 
under the addition of mass to maintain stability (Chandrasekhar 1939). Slow, spherically 
symmetric AGB winds may therefore be the reaction product of core contraction, whereon 
they would be the three dimensional analog of a rocket's exhaust. Similar winds occur 
during one stage of star formation (e.g. Shu & Lada 1990), while a protostar is separating 
itself from its molecular shroud. At a more advanced stage of the process proto-stellar 
winds usually seem to develop into bipolar flows, and the resulting T Tauri stars often 
have an annular disk of rotating matter. Bipolar flows occur in the last stages of AGB 
mass-loss too, and characterize most if not all proto planetary nebulae (Hrivnak 1996). 
Thus spherically symmetric winds seem to be a feature of active gravitational contraction 
when angular momentum constraints are not important. 

The simplest systems with significant mass loss are those with little complexity added by 
rotation or magnetic fields, such as the AGB stars. We therefore explore the basic 
hypothesis here that the momentum flux required to drive a spherically symmetric AGB 
wind can be attributed to the excess to-virial-stability momentum flux generated by 
gravitational contraction after it has been transmitted through a star's center to become an 
outwardly directed flux. A similar inversion in direction occurs with standing waves in a 
spherical cavity (Landau & Lifshitz 1979). The central question posed by this paper 
concerns the proper treatment to accord this momentum flux, which is usually assumed to 
be cancelled by symmetry, without considering how that is achieved in detail. 

Section 2 exhibits the present theory in outline form by calculating the momentum flux 
generated in the quasi-static contraction of a homogeneous sphere. A general formula for 
the flux from a contracting polytrope under the addition of mass is derived in section 3, 
and applied to AGB stars in section 4. A simplified theory is applied to the Sun in section 
5 to explore the relationship of this approach to the momentum flux in the solar wind. 
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2. the secular contraction of a homogeneous sphere 

The electron-degenerate core of an AGB star secularly increases in mass as helium 
"ashes" from shell hydrogen burning are added (e.g. Iben 1987): the core contracts to 
maintain virial equilibrium. Our objective is to quantify the resulting inwardly directed 
momentum flux for comparison with what is needed to sustain mass-loss. Both for 
simplicity, and to exhibit the essential features of the scenario quickly, our program is 
carried through in this section by treating the core as a homogeneous sphere. 

Consider first a spherical shell of radius r and mass dm = 4 nr^pdr, with every particle 
in circular orbit about a central mass, M{r), so no particle has a radially directed velocity 
component. The change in the potential energy of the shell, At2, as r^r-Ar in an 
infinitesimal contraction is 



= GM 



r 1 



yr— Ar r J 



dm = dmv^gj{r) 



^ Ar ^ 
Vr-Ary 



SO its kinetic energy immediately after a radial contraction Ar is to first order 



Tsheii = ymv^otif) + dmvL{r) 



' Ar ^ 
Vr-Ary 



= ymvUr)[l + 2Ar/r]. 



The radially directed change in the velocity of every constituent particle is therefore 



where is the fractional change in its radius. The resulting radial momentum flow has a 
magnitude A P = dm Av (just as that resulting from the free-fall of a shell from infinity is 
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AP{dm, r <- oo) = dmv^J, since a net inward movement of mass constitutes a local 
momentum flow, despite the fact that spherically symmetric contraction does not change 
the vectorial sum of momentum in the system. This use of A P emulates the ubiquitous 
L/c formulation (Salpeter 1974; Knapp et al 1982) for quantifying the total momentum 
flux carried by a star's luminosity in every radial direction. 

We now extend this analysis to the simultaneous contraction of every shell of a stellar 
core in hydrostatic equilibrium. The resulting inwardly directed momentum change of 
each constituent shell is part of a coherent radially inward momentum flow. This cannot 
be cancelled by interactions within a mass distribution in hydrostatic equilibrium, where 
every other momentum component is by definition balanced, as Parker (1958) long ago 
emphasized. And inwardly directed flows have no possibility of being exhausted by work 
against gravity. The momentum flow must therefore reach the center. 

The inwardly directed momentum increment generated by a proportional contraction of 
every radius of a homogeneous sphere by p is thus 



3^^ A 



ATI r^'o ^ \A 2 , of^s ATtr'pG 



111 



47tr'pdr = ^(5Mv^^^ (1). 



3V2 



Here Rg is the surface radius, M the total mass, and v^^^ = ^ 2 G M / the escape 
velocity. When contraction is part of a continuous, quasistatic process, differentiation of 
AP with respect to time gives the momentum flux 



dP 3 V 2 M dRc 

— = V (2) 

dt 8 Rs dt ^ ' 



through the center engendered by a contraction rate dRg/ dt 
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To progress further we need an expression for dRg/ dt. Our AGB core is represented at 
this point by a homogeneous sphere with the same surface radius as the core, whose 
luminosity, Z*, is related to its mass, M^, in solar units by the relation 



LJ L„ = 52000 (M, / 0.456) 



(3) 



deriving from stellar structure calculations (Boothroyd & Sackman 1992). This permits 
to be obtained from L* by dividing it by the energy released per gm of hydrogen converted 
to helium, E^, which after allowing for the mass fraction in hydrogen, X, gives 



dMJdt = Lj{XEu) 



(4) 



(Kwok 1987). The core contraction rate consequent on the mass increment rate is got 
by differentiating the relation between & given in Hansen & Kawaler (1994; pi 27) 
as a fit to the results of exact calculations for degenerate white dwarfs 



Rs/R, = 2.02[l -{M J M^y'Y" {M J M^y 



1/3 



(5), 



where i?i = 1.117*10 ^ ii^^ R^, M^ = 1.456 (2 /^J^M„, and /j.^ the mean molecular weight 
per electron, so 



dRs 
dt 



2. 02 R, dM^ 



3M„ dt 



{i-(M,/M.r} 



1/2 



(6). 



The minus sign signifies that the polytrope contracts on the addition of mass. Substituting 
(4) & (6) into equation 2 gives the desired expression 



The momentum flux needed for lifting mass from the photosphere to the radius at which 
dust condenses is ~mv^^^. For the canonical Mira of Humphreys et al. (1996), with an 
effective surface temperature T^^ = 3000 K, period of 332 days, Rg = 244 Rq, expansion 
velocity v^^p ~ 10 km s""", and m = 1.8 10"^ Mq yr""", L, = 4340 Lg, which from equation (3) 
implies an = 0.539 Mq. This has a photospheric v^^^ = 39.5 km s ^ so for M = 1.0 Mq 
(core + envelope) 

mv,,, = 0.081 (Z./c) (8). 

This should be compared with the contraction flux calculated from equations 4 & 7 with 
X= 0.7 and £^ = 6 10^^ ergs gm-i, of 

— = 0.850 (L,/c) (9). 
dt 

Clearly the momentum flux engendered by core contraction is in this case an order of 
magnitude larger than that directly needed to sustain the modest mass-loss rate in (8). In 
the absence of an efficient central momentum cancellation mechanism, this flux passes 
through the center to become in turn an outwardly directed flux that is transmitted from the 
core to the stellar envelope, to promote and sustain its characteristic pulsation mode. This 
in turn leads to a shock distended stellar atmosphere, and thus to photospheric mass loss 
(Bowen 1988). 
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3. Contraction Momentum from Polytropes 

White dwarfs have a polytropic mass distribution (Chandrasekhar 1939). Our next 
objective is accordingly to quantify the momentum flux generated by a polytrope under 
changes to its radius, Rs, and mass.M. We begin with the gravitational potential energy, 
Q., oi a stable polytrope of index n (Chandrasekhar, p101), which is related to its total 
kinetic (or thermal) energy content, T, through the Virial Theorem such that 2T + Q. = , 



5 — n 



GM 
J Rs 



SO 



T = 



in = i 

2 2 



5 — n 



= -MKk.c = 2 """" 



with A „ = 



r 3 A 



1/2 



The rms velocity, v^, of particles (for the moment assumed to 



all be of one kind) is related to their mean velocity, v, in a polytrope by 



v = Av = A, A, V 

'^rms ^ ^ n * esc 



(10), 



where X = 0.921 when they are from a Maxwell-Boltzmann distribution. Let us define P, 
the total "directionless" or scalar momentum content of a polytrope in virial equilibrium, as 



p = Mv = MXXv^ 



(11), 



which relates to the usual expression between momentum and kinetic energy via 
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The value of A depends on the motion of constituent matter. For an n = polytrope (the 
homogeneous sphere) with every constituent particle in circular orbit about the center, 
following the development of (1) above, v is 



V = 



1 7 



M 



? r ^. 2 ^ 1 "((GM{r)\'\ , ^ 

^rot v) ^ J^rp dr= — ^— 4 TT ry dr = 

i M i\ r J 



3 V2 



(13), 



and since (10) defines v = A v^^^ , A is then 



= 0. 559 



The homogeneous sphere has utility in treating the momentum content of polytropes, as it 
is easily integrated. The results are then generalizable to an arbitrary polytrope of index 
n, since this can always be mapped to the equivalent homogeneous sphere with the 
same total binding energy and mass, by evaluating its surface radius, ^, via 



T = 



A o A 



5 - n 



GM^ 1 



GM' 



5J R. 



S,o 



SO 



i?c — — 

S,o c S^f^ 



R 



S,n 



(14). 



It follows that A„ v^^^ „ = A„ v^^^^, whereon the momentum content of the polytrope, P„, is 



P„ = Mv„ = MAA„v,,,,„ 



Mv„ = P„ (15) 



This gives the exceedingly useful result that the P of a polytrope of index n is readily 
derived from that of its equivalent homogeneous sphere. 



The second use of P comes from its differentiability: that with respect to radius gives 
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dP = AAq M 



dv. 



dRs = - 



1 3 V2 M 



1 PdR^ 

2 R, 



'^=-^PP (16), 



dRs 



2 S Rs 



while a similar differentiation of the 7 of a poytrope gives 



dT = 





(17). 



Both dP and dT are positive when account is taken of the decrease in radius attendant 
on contraction. So (16) & (17) show that dP = 0.5{P / T)dT is proportional to the 
change in the kinetic energy of the polytrope. It should be noted that while (16) gives the 
change in the momentum content of a polytrope in virial equilibrium, dP = -0.5 p P, {^) 
gives the change dP = - P P resulting from an instantaneous contraction before any 
contraction energy is exported. The factor of two difference in coefficients 



is simply the ratio between the change in retained kinetic energy of a polytrope in 
equilibrium and the change in its potential energy under contraction, as half of the 
contraction energy in non-relativistic regimes must be exported to maintain virial 
equilibrium. In consequence a summation over the inwardly directed momentum 
generated by quasistatic contraction through -AR^ has a magnitude of twice the change 
in the equilibrium value of P for the polytrope: this is exactly quantified by the 
differentiation of P with respect to R^ 



O.sp P _ I _ AT 



PP 2 Ml 



The rate of change in the P of a polytrope under changes to both its mass and radius is 
obtained by differentiating (11) 
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dP dp clM dP dR. 3v dM v M dR. 

— = + ^ = ^ (18). 

dt dM dt dRs dt 2 dt 2 R^ dt 

The sign of the second term in (18) is cancelled during contraction by that of dR^I dt, so 
the terms add and will be seen below to be of comparable size. But the AP generated by 
the addition of AM is at most AP = 2 (AM I M)P (Appendix A, equation A2) when AM is 
brought from infinity, whereas (18) shows the change to the Pof the polytrope alone 
under the addition of mass at constant R^ is (3/ 2) (AM I M)P: a doubling of the first 
term in (18) does not in this case provide its contribution to the total momentum generated 
by contraction. Thus the mode of mass addition needs more scrutiny. 



Mass is added to an AGB core from its periphery as it undergoes quasistatic contraction, 
so it is appropriate to consider the AM mass addition as coming from a thin shell at a 
radius (P^+Ar) with a prior kinetic energy AT = 0. 5 AMv,^„,(i?5+ Ar) and a prior 
AP(AM) = AMv,„,(P5 + Ar). This can be rewritten in the case of an n=0 polytrope as 



AP(AM) = —AM 
^ ' 2 



' 2GM{r) ^"' 
V Rs + ^r J 



4 (^M_ 
sl M 



(19), 



where P„ = 



3 V2 



M V,,, from (1 1) & (13). The difference 



5P = 



fAM^ 




rAM' 








1^ = - 




[ M . 




I M , 


> 6 


I M ) 



(20). 



between the net change AP = 1.5( AM/M)P of the mass-augmented polytrope in 
virial equilibrium at constant R^ after the addition of AM from (18), and the AP of (19) 
which is the momentum content of AM while it is still at (Ps+ Ar), is the increase In the 
polytrope's P due to homologous contraction as m[Rs) M{Rs-Ar) and 
AM (P^ + Ar) AM{Rs). Since (20) originates solely from contraction, it too represents 
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a radially directed momentum flow. It follows from (18) & (20) that the scalar sum of 
"freshly generated" momentum from a polytrope in virial equilibrium, consequent on 
changes to its mass and radius, is 



AP = -V AM - - V 
6 2 



Ar 



(21). 



Differentiating (21) to cast it as a flux, and doubling its coefficients to get the inwardly 
directed momentum flux, gives the basic momentum flux relation for any polytrope 



dP I _dM _ 

— = — V V 

dt 3 dt 



dRs 
dt 



(22). 



Equations 4 & 6 are then used in conjunction with (22) to predict the momentum flux 
arising from the quasistatic contraction of an AGB core. 



4. Application to AGB stars 

This section deals first with the momentum flux from the core of an AGB star, then with 
what is needed for mass-loss from its photosphere, before finishing with the results 
accruing from applying this approach to understanding the winds of AGB stars. 

4.1 : on the momentum flux from core contraction 

The whole of a degenerate core contracts in reaching virial equilibrium on the addition of 
mass, and half of the resulting contraction energy is exported (less is exported from highly 
relativistic cores). And though the instantaneous scalar momentum change AP oc AO., 
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the loss of -50% of the contraction energy does not of itself imply an exactly proportionate 
export or cancellation of 50% of AP, since this outcome is mediated by (i) equipartition of 
energy between particle species, and (ii) the core is supported against collapse by 
electron degeneracy rather than thermal pressure. We look at these factors next. 



(i) equipartition of energy 

Suppose for simplicity a non-degenerate core composed entirely of carbon-12 nuclei, so 
there are 6 electrons for every nucleus. Under equipartition of energy the kinetic energy 
of all nuclei in a polytrope in thermal equilibrium is 



T 
1 



I M_ 

7 y 



n esc 



m, 



Mv' 



whereon 



1 + 6 

^7 J 



SO for the nuclei alone P_ = 



f 




f 








^mi2+ 6 




V 



1 m, 



1/2 



7 m,T + 6 m„ 



P (23), 



while the contribution from electrons is = 



6 6 OT„ 



y7 mi2 + 6m, y 



The ratio of the scalar momentum content carried by electrons to that by nuclei is 
6(m,/ mi2)"^ = 0.04, and the ratio P„„^ / P = Vl/7 = 0.38, or {/i I jijY'^ whenever the 
number of electrons per nucleus is ~i^j/2, where ^ and ii, are respectively the mean 
molecular weight and the mean molecular weight of ions. Clearly the equipartition of 
energy between particle species results in the scalar momentum content of a polytrope in 
equilibrium being much smaller than if it were composed of just one species. 
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(ii) the scalar momentum content of nuclei in a degenerate core 

Since core nuclei are supported against gravity by electron degeneracy rather than 
thermal pressure, their kinetic temperature is markedly lower than that needed for virial 
equilibrium under thermal pressure. And their temperature, r„„^, is necessarily less than 
the = 10^ K that triggers conversion of He^ to C^^ during shell-hydrogen burning 
phases of an AGB star. The cores are therefore nearly isothermal. The mean kinetic 
energy of a carbon-12 nucleus is (1 / 2) = (3 / 2) A: r„„^, so the scalar 

momentum content, P„„, = N v„„, = Nmu^ v^^„^ of an entire core is 



where k is Boltzmann's constant, A^^ Avogadro's number, and N = M N^/ /Hi the total 
number of ions in the core. Since nuclei have a Maxwell-Boltzmann velocity distribution, 
A is in this instance 0.921. 

(Hi) the scalar momentum content of the electrons 

We derive an approximate value for the scalar momentum content carried by the electrons 
of a degenerate core, when this is treated as a homogeneous sphere next: this follows the 
treatment Hansen & Kawaler (1994, p124ff) accord the Chandrasekhar mass limit. The 
electron number density, n^, of a completely degenerate, zero-temperature, non- 
relativistic sphere whose Fermi-momentum = x^(m^c), is 



where h is Planck's constant. But n^ = p N ^1 ii^,so for a homogeneous sphere is 



P 



(24), 
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9M 

32 TT^I^ 



1/3 



J_ 



(26). 



The scalar momentum content of the electron distribution, after substituting (26), is 

= 1.^3 r.. 8. ^ Sn^(h_^^9MN^Y' (27). 

This provides an upper limit to P,, since all other polytropes have a larger central density, 
which increases the mean electron energy in the center and reduces their contribution to 
. When (27) is evaluated for our prototype Mira and P„„^ is calculated on the same 
basis from (24), P^ I P^^ = 0.08. Though most of P resides in the ions, this changes as 
the core mass increases, as a result of the concomitant decrease of Rg on (27). 



The scalar momentum content of a body in virial equilibrium is the sum of contributions 
from its ions and electrons, P,„re = Pnm + Pe ■ This is smaller than the P defined by (11) 
when the existence of more than one kind of constituent particle is ignored. For our 
prototypical Mira core of §2, as it contracts through - A/? under the addition of mass 

AP,„,, = [{AP„^I AP) + {APJ AP)}AP - {0.1012 + 0.0081} AP - 0.11 AP, 

where AP is the total change in P engendered by core contraction from (22). It is clear 
from this calculation that -90% of the momentum generated in quasistatic contraction is 
either exported from or cancelled in the core in maintaining equilibrium. 



The momentum flux equation (22), on substitution of (4) & (6) becomes 
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dt 



, 2. 02 R, M, 
1 + ' 



R 



-4/3 



1 + (M, / M„) 
{l-(M./M.)-}'" 



(28), 



where is the mass and R^^ the surface radius of the core. This is evaluated with 
V = A A„v,,,_„ = (3a/2 / 8) ^5 / (5 - n) v,,,_„, in accord with (10) & (15). It follows 
from the discussion above that a fraction e = 1-{P„^^ + P,) I (2 P) of this radially 
directed momentum flux is surplus to the stability of the core. 



Table I lists some representative values of the momentum flux generated by different core 
masses when (28) is evaluated with X= 0.7, n = 3/2, and r„„^ = 10^ K, together with (3). 
The flux increases much more quickly than the luminosity with M^, as is evident in its 
order of magnitude increase down the table when expressed in units of L, / c of the 
stellar luminosity divided by the velocity of light, which already discounts the associated 
increase in luminosity. This portion of the flux increase is due in part to a modest increase 
in the fraction available for export as e increases with M^. But most of it is due to the 
accompanying decrease in R^^, which causes a rapid increase in the second (radial 
contraction) term in the main bracket of (28), whose relative value increases by a factor of 
5 down the table. In real cores this effect is enhanced further by n increasing with M^. 



4.2: momentum flux requirements of a slow wind 



The second part of the stellar wind problem is to determine what momentum flux, P , is 
implied by the occurrence of mass-loss from a photosphere in the ideal circumstance 
where none is cancelled by shocks in the distended atmosphere: more is of course 
needed to compensate for these losses. The smallest momentum input enabling mass 
loss from a photosphere corresponds to mass leaving at a velocity sufficient to reach the 
radius of dust formation, r^. When it is further assumed that = oo, we have the 
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impulsive limit P = m v,,,^, where v^,^ ^ is the escape velocity from the photosphere. But 
Bowen's mass-loss models have mass leaving the photosphere at a smaller radial 
velocity in a more massive mass flux, which requires a larger P. Quasistatic contraction 
has the opposite bias to that of the impulsive limit, as it maximizes the quantity of radially 
directed momentum generated during contraction. We therefore use its formalism next to 
place an upper limit on the P required to sustain mass-loss from a photosphere. 

During quasistatic contraction of a polytrope from to Rs, the change to its Pis 



P[Rs) - P[r,) = 



3 V2 TA, ^ 



M sjlGM 



d J 



(29) 



from (11). Similarly, when the polytrope has instead a mass M - dM 



P{Rs) - P{r,) = 



8 [^o^ 



[M-dM) J2G{M-dM) 



d J 



(30), 



so the difference between (29) & (30), retaining first order terms, is 



A/'(AM) = |^[p(i?,)-P(r,)] (31). 



Finally by adding to (31) the term from (20), which allows for the momentum-debt of 
subtracting mass from a polytrope while Rs = r^, we get the quasistatic change to the P 
of a polytrope when it expands from R^ to with mass M, loses dM, and then contracts 
to Rs- These operations on the equilibrium values of P for a polytrope represent just half 
of the total scalar momentum change implied by this quasistatic cycle, whose net result is 
the loss of dM at zero velocity from r^. The momentum flux required to quasistatically 
power a wind (the time-reversed inverse of contraction) is therefore 
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AP(AM) = \^[9P{Rs) - 8P(r,)] 



(32). 



In the interesting limit of ^ oo , (32) expressed as a flux reduces to 



P = 3 



3 V2 TA, ^ 



m v,,,,p (33), 



which gives P = 1. 9016 m v,,^ ^ for an n = 3/2 polytrope. 



The two momentum flux bounds on mass-loss from the photosphere as r^^oo are 



mvescp < P < 1-9 mv,,^p 



(34). 



These are both made smaller by taking detailed note of r^, as mass is then lifted through 
a smaller height against gravity before radiation pressure on dust takes over, whereon 



mv 



esc,p 



1_ i^s 



V 



< P < V2m 



esc,p 



8 ir. 



(35). 



When a model dependent is related to a usual dust condensation temperature -1000 

K in an oxygen-rich shell (e.g. Humphreys ef a/1996) via = Rs[Teff 1 1000)^'^ and 
» 3000 K, is 15.6 Rs- This gives our 1 Mq Mira of §2 the momentum flux bounds 



0.747mv,,,,^ < P < 1.420mv,,,,p (36). 



The P needed for lifting mass from Rg to is set by these radii. These depend on T, 



eff 
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whose value is correlated with both m and progenitor mass. Most Miras have almost 
translucent shells with T^^ ~ 3000 K, while heavily shrouded objects have much lower 
temperatures, though the actual T^^ are moot. Le Sidaner & Le Bertre (1996) assign T^^ 
of 1800 - 2200 K to OH/IR stars, whereas a very similar shell-modelling study by Lepine, 
Ortiz & Epchtein (1995) assigns many of them a T^^ ~ 1250 K. WX Psc is an intermediate 
case. It is a heavily shell-reddened star with an M9.5-10 spectral type (Lockwood 1985), 
and so perhaps may be accorded a ~ 1600 K (Dyck et al 1974). We expect the 
reddest OH/IR stars to have the lowest T^^ , which is supported by the correlation between 
spectral type and expansion velocity (Lewis 1991). In practise the stars with the largest m 
have the coldest photospheres, which means they raise their mass flux through a 
proportionately smaller radial range to reach r^. These objects therefore have the innate 
capacity to support the largest m when mass-loss is a momentum bounded process. 
Table 2 illustrates these trends for a 1 Mq star. Its last column shows that a 3.4 times 
larger momentum flux is needed to drive the same m from a star with T^^ = 3000 K as 
from one with T^^ = 1 500 K. 

4.3: results for AGB stars 

We saw in §2 that contraction readily provides for the needs of an m = 1 .8 10 Mq yr 
loss from our prototype Mira with T^^ ~ 3000 K. But an important fraction of any 
momentum flux is cancelled in photospheric shocks, and some may also be cancelled at 
the center, so a considerable extra capacity for generating a flux, over and above the 
direct needs of the stellar wind, is a mandatory feature of any theory. We work hereafter 
from the basis that the minimal flux is specified by the quasistatic upper bound of (35). 
The most demanding oxygen-rich AGB stars to be accounted for have an m ~ 10 Mq yr 
(Knapp & Morris 1985). For an assumed = 1500 K, the quasistatic bound for 
sustaining this wind from a 1 Mq star with L = 30,000 Lq is 1.76 L, / c; if T^^ = 1750 K , 
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the bound is 2.47 L, / c. These bounds are both suitably smaller than the P = 5.51 
L, / c listed in column 5 of Table 1 as being generated by core contraction at this 
luminosity. Since the AGB stars with the largest m may have L = 50,000 Lq, and so may 
generate a commensurately larger P, we conclude that core contraction is easily able to 
provide for the momentum flux needs of every AGB wind. 

An immediate result of attributing the AGB wind's momentum flux to gravitational 
contraction is its consequent sensitivity to M^, and thus to L. AGB stars with small 
progenitor masses, and so relatively small luminosities, absolutely can not achieve the 
largest observed mass-loss rates. This has already been forced on our notice by 
observations of circum-stellar masers that show a marked decrease in the galactic latitude 
distribution of OH/IR stars with increasing shell opacity, and hence m, when these 
parameters are assessed from their IR colors (Lewis 1987). An update of this evidence is 
presented in Figure 1 for the IRAS sources with detected OH masers: Figure 1 shows an 
increasingly marked paucity of high latitude objects when (25-12) |im > -0.4, as shells get 
increasingly redder and more opaque. Column 6 in Table 1 illustrates the same point 
numerically in documenting the rapid increase in the maximum m with M^. This is 
obtained by equating the P calculated with (28) & (3) to three times the quasistatic shell 
bound computed for T^^ = 1 500 K: the factor of 3 here is equivalent to the assumption 
that just one third of the generated momentum flux survives cancellation in atmospheric 
Shockwaves to support mass-loss at all progenitor masses. On this basis an m » 10 Mq 
yr implies an L » 10,000 U and ~ 0.65 M q. 

Some predictions flow from the applicability of our premise. Thus a sine qua non for the 
occurrence of mass-loss from an AGB star is that its core must be contracting. Yet during 
a thermal pulse the outer helium shell is stripped from a core as it burns, causing a net 
loss of core-mass. The core must then expand to maintain virial stability through the 
event, and so becomes a net consumer of momentum, which implies that m ~ during 
this phase. Indeed CO observations of U Ant, S Set, & TT Cyg show substantial 
decreases in their mass-loss rates in the wake of a thermal pulse (Olafsson et al 1993). 
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Part of this is attributable to temporarily smaller core masses, but most is due to the drop 
in luminosity that follows on a thermal pulse (Wood & Zarro 1981), which implies a 
correspondingly lower rate of mass accretion. On the other hand post-AGB stars at the 
termination of mass-loss, following the dissipation of most of the stellar envelope, 
continue to burn hydrogen as rapidly as ever, and so continue to increase their core- 
masses. Mass-loss should therefore continue in some form as the most efficient mode for 
the dissipation of P, even if the nature of the wind and its acceleration changes. "Fast 
winds" have long been observed to emanate from the central stars of planetary nebulae. 

There are natural extensions of these ideas to stars with angular momentum (and/or 
strong magnetic fields), since the second sine qua non of the mechanism is point 
symmetry rather than spherical symmetry in accretion. A degenerate core spins faster as 
it grows in mass and contracts, so the angular momentum of material will at some stage 
affect the local accretion rate. This is necessarily slower at the equator than the poles, so 
an erstwhile spherical symmetry in the production of a momentum flux develops axial 
symmetry. This sculpts the wind, which refashions its environs. Material is then swept 
away from the poles more quickly than from the equator as there is a larger flux there, 
which promotes the probability of both equatorial disks and bipolar flows. These factors 
operate in every AGB star (and in every mass accreting object), not just those with stellar 
companions. The bipolar structure of most planetary nebulae (Zuckerman & Aller 1986) 
and proto planetary nebulae (Hrivnak 1996) is therefore an expected result of the 
applicability of this theory. It may also be applicable to star forming regions, since these 
mirror the same progression from spherically symmetric winds to bipolar flows. 



5. The Solar Wind 

Red supergiants have slow dusty winds rivaling those from AGB stars without having 
degenerate cores. By analogy they must also generate momentum fluxes, which in this 
case occur as much as a byproduct of stellar fusion as of contraction. Since the central 
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zone of a star contracts to maintain tliermal pressure as its particles are consumed in 
nuclear reactions, this effect may suffice to generate the P required to sustain mass loss 
from red supergiants, a supposition that needs testing with a stellar code. We use the 
idea here with a simple polytropic model of the Sun to see how it relates to the solar wind. 

Let us model the Sun with an n = 3.25 polytrope. This maps into the homogeneous 
sphere with the same binding energy via (14), so we can predict a contraction rate by 
assuming it maintains a constant particle number density as fusion reduces its particle 
count. The assembly of each He^ nucleus reduces this by 5, so dNI dt = VaN^ {l^ I 
The rate of volume change of the homogeneous sphere is then 



d{%KRl) _ 5 



dt 



( T \ 



where ^ is the mean molecular weight and p = M„ I r], whereon 
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(37). 



This relation is unchanged by mapping between the sphere and polytrope. Substitution 
into (2) gives the momentum flux, implied by homologous contraction of the solar 
polytrope, of which a fraction e is available for export following an argument like §4.1 (I), so 



dP 
dt 
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(38). 



The numerical value comes from using }i = 0.613, and 8 = {1 - (X/2 + V73 + Z/3) / 2} ~ 
0.774, in which X , V, & Z, the mass fractions of H, He and heavy elements, are 
respectively 0.7094, 0.271 & 0.0196. For v^^^ j ^j = 617 km s P„ is 7.0 10^1 gm cm s-2 
(i.e. 0.055 L, / c). The P powering a quiescent solar wind at a velocity -400 km s"^ and 
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an m « 10 yr -1 (Cassinelli & MacGregor 1986) is from (35) with n = 3/2 

PoMnd = ^2m{^n I ^o) [v...,3.25 + 400] 9/8 = 1.219 10'*^ gm cm 5"^ (39). 
Thus P„ is -57 times larger than the quiescent solar wind needs. 

The momentum flux required to sustain the totality of components subsumed into the solar 
wind is uncertain, though it is likely to be 2-3 times larger than (39). Thus apart from the 
flux originating from flares and other impulsive outbursts, Thomson ef a/(1995) discuss 
evidence linking an additional impulsive component of the wind to standing waves in the 
solar structure. These are omnipresent, encompass the whole surface, and may on their 
own more than double the particulate flux emerging from the Sun. Even so the P„ from 
(38) is still an order of magnitude larger than the whole of the solar wind needs. 

The present calculation is simplistic since the radius of the Sun increases slowly as it 
evolves. This reflects the influence of other factors on its structure, including a very 
modest rise in central temperature, which decreases the available momentum flux and 
opposes the tendency for the whole mass distribution to contract. But the basic structure 
of the Sun is determined by hydrostatic equilibrium, when the tiny perturbation caused by 
fusion is ignored, so the contraction momentum flux must either be transmitted to the 
surface and/or used to expand some of the layers overlying the nuclear burning zone. 
The differential importance of these options needs elucidation from a stellar code. 

It is instructive to compare P„ with the magnitude of the momentum cancellation rate 
inherent in stellar fusion. The cancellation rate from the dominant, non-resonant, p-p 
reaction can be approximated by assuming it occurs at the Gamow peak energy, 
e„ = 1. 22 (zf Z2 ;Ur T'e)^' where /n^ is the reduced mass of the interacting particles (1,2) 
in amu, in units 10^ K, and e„ in kev. The momentum cancelled per He atom is 
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AP(He'*) = 2 ^m^e^ , which leads to a solar p-p cancellation rate when = 15 of 



dP N. I 

dt 2 ^ ^ " 



= 2.43 10^^ gm cm s'' (40). 



Since (40) is much larger than (38), there will normally be NO equality between the flux 
generated by contraction and that cancelled in fusion. In actuality momentum cancelled 
in fusion reactions is returned to the stellar structure as fusion energy is thermalised. 
Since fusion is much the most efficient stellar channel for cancelling momentum, this 
result suggests the contention that the scalar momentum content of a stellar structure is 
statistically constant so long as its mass distribution and its ionised mass-fraction remain 
uncfianged. 



6. discussion 



Gravitational contraction always generates an inwardly directed momentum flux. The 
magnitude of this flux is quantified here for the quasistatic contraction of a polytrope, 
which is well suited to describing the secular contraction of an electron degenerate core 
under the addition of helium ashes: it is also applied to the contraction of the nuclear 
burning zones of the Sun as its particle number count changes with fusion. The 
momentum flux generated by gravitational contraction is, on the assumption that it is all 
transmitted through the center to become in turn an outwardly directed flux, more than 
sufficient to provide for the momentum requirements of slow, m ~ ^0'^ Mq yr AGB winds 
on the one hand, as well as those needed by the m = 10 ''^ Mq yr"^ solar wind on the 
other. A mass efflux is probably a general concomitant to gravitational contraction 
whenever a system has insufficient resources for cancelling its radially coherent 
momentum flux. An alternate resource is pulsation, which allows radially directed 
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momentum to be thermalised in atmospheric Shockwaves. The known correlation 
between pulsation and mass-loss is therefore explicable as a correlation between 
independent pointers to the existence of a significant, radially directed momentum flux. 

The application of this causality to the specific case of AGB stars explains the extreme 
sensitivity of their maximum mass-loss rates to progenitor mass. Indeed the momentum 
flux from their degenerate cores has a magnitude 1-10 times that of the radiant flux in 
Table 1 , which necessarily has implications for the structure of their envelopes. This is 
certainly one of the reasons Red Giants are giants, as their contracting cores are sources 
of a significant momentum flux that has to be absorbed and/or dissipated by the envelope. 
Since the essential feature underwriting the ability of the contraction mechanism to 
generate a momentum flux is point symmetry about the center, any important axially 
symmetric feature of the core, such as its rotation or magnetic field, can influence the local 
accretion-rate onto the core. This in turn impresses an axial signature on its contraction 
momentum flux, which will lead in turn to large scale structural features, such as bipolar 
flows and equatorial disks. The prevalence of bipolar flows in both proto planetary 
nebulae and planetary nebulae is thus an explicable consequence of the export of the 
contraction momentum flux from a rapidly spinning core. 

These conclusions follow from the realization that there are several ways in which a 
spherically symmetric momentum flux can be treated, so that the overiding principle of 
conservation of momentum is satisfied. One mode is for it to be explicitly cancelled at the 
center of symmetry, such as happens with single bubble sonoluminesence. But no 
analogous mechanisms are available to stellar systems. A second mode is for it to be 
thermally cancelled near the center, which is difficult within the constraints of hydrostatic 
and thermal equilibrium obtaining in stellar structures. It is, however, reasonably easy to 
achieve some cancellation by disturbing the local equilibrium at the periphery of a star. 
This can be facilitated by pulsation, as happens with the ejection of energy into the solar 
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wind by standing waves in tine Sun (Tliomson et al 1995), and is more generally 
evidenced in the very high temperatures of stellar coronae. Mass-loss in a wind provides 
a third mode for disposing of a radial momentum flux, that inherently maintains symmetry. 

Our conclusions are premised on the supposition that the majority of the contraction 
momentum flux generated in the core is exported, whereas it is certainly reasonable to 
expect that in the absence of any other mode of disposing of it, it is thermally cancelled. 
This is the usual if little explored expectation. So the central question posed by the 
present exercise is the proper treatment to accord the contraction momentum flux. 

Since the flux originates from a mechanical response of the stellar structure, it implicitly 
has a dynamical time scale. Its complete cancellation at the center, however, would 
associate it in most stellar contexts with the much longer photon diffusion timescale: these 
differing timescales can only be reconciled by the adjustment of the equilibrium thermal 
gradient around the center to the needs of cancellation, together with a perpetual 
mechanical oscillation of the whole structure. Yet this class of solution immediately 
requires the momentum flux to flow repeatedly through the center. Unless the momentum 
flux is immediately and completely cancelled at the center on its first arrival there, and 
there are no efficient mechanisms in any star for accomplishing this, its passage through 
the center to become in turn an outwardly directed radial flux mandates its transmission to 
the rest of the stellar structure, so that every other mode of disposing of a momentum flow 
needs exploration. The most obvious are (i) the consumption of the momentum flux in 
work against gravity, the inverse of the process from which it originated, and (ii) its export 
from the entire structure as a wind. 

It is worth pondering further on the situation posed by an AGB star when contraction is not 
deemed to be the source of the momentum flux driving its wind. In that case the wind's 



29 

flux must originate from some shell concentric with the center that lies beyond the core. 
And since action and reaction are equal and opposite, the launching of any momentum 
towards the outside implies the simultaneous and equal launching of a momentum flow 
towards the center. But the basic structure of a star is set by hydrostatic equilibrium, the 
more so as its center is approached, so the inward flow can neither be cancelled, nor 
absorbed by doing work against gravity, and so should reach the center. For consistency 
in this context, the inward flow has to be completely cancelled by the mechanism held to 
dispose of the central contraction momentum flux. Acceptance of this scenario therefore 
requires (i) an efficient, as yet unrecognized, central momentum cancellation mechanism; 
(ii) a mechanism for generating a radially directed momentum flow from a larger radial 
zone in the star; (iii) additional mechanisms for explaining large scale features such as 
bipolar flows (in AGB stars appeal is currently made to interaction with a close companion 
star, despite their rarity); (iv) auxiliary mechanisms for explaining winds from stars and 
proto-stars; (v) auxiliary hypotheses for explaining the large scale structural morphology 
of other contraction contexts. The need for an array of mechanisms is resolved by one 
application of Occam's razor. It is surely more reasonable to expect most contexts in 
which mass loss occurs in propinquity to gravitational contraction to be influenced by the 
inescapable generation of a contraction momentum flux there, than to invoke an array of 
mechanisms to replace it. My particular conclusion is that the ultimate source of the 
momentum flux powering every slow AGB wind is the momentum flux generated by its 
contracting core. 

An early crude version of this paper was prepared while I was on leave at the Max Planck 
Institute fur Radioastronomie; I thank Professor Wielebinski and Dr. Tom Wilson for their 
hospitality and a relaxed environment. This work is supported by the National Astronomy 
and Ionosphere Center, which is operated by Cornell University under a management 
agreement with the National Science Foundation. 
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Appendix A 

If the whole of the potential energy change of a homogeneous sphere under the addition 
of a mass, AM, is available to increase its kinetic energy at the moment of addition, then 
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M )' 



When p ^ p + Ap under this change, the kinetic energy of a spherical shell of mass 
dm= Aitr^pdr and radius r in circular orbit about a central mass, M(r), becomes 

T^ueii = ^dmvl{r) + dm[2\^\AM{r) I M{r)] = ^ dmvi{r)[l + 4 {AM {r) / M {r))] , 



SO the change in the velocity of dm as a result of adding AM(r) can be thought of as 

Av(r) = 2v,„,(r)(AM(r)/M(r)) (A1). 

The resultant change in P, the scalar sum of the momentum content of the 
homogeneous sphere defined by equation 11, when AM{r) = 4 n r^Ap/3, is 
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so from (11) & (13), AP=2P{AM/M) (A2) 
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Table 1 : change in the contraction momentum flux from a 3/2 polytrope with luminosity 
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5.506 
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12.898 
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t this upper limit to m is obtained by equating P to three times quasistatic limit of (35), 



with T,ff = 1500 K 
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Table 2: 

sensitivity of tlie momentum flux needs of photospheric mass-loss to r, 
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39.5 
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2500 


351.9 
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3.14 
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23.0 
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2.45 
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2.76 
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t the P/ m ratio is calculated from the quasistatic bound of equation 35 
for a 1 Mo star with L = 3440 Lq and n = 3/2 . 
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Legends: 

Figure 1 : The galactic latitude distribution of IRAS sources with detected OH masers 
plotted against (25-12) |xm IR color index. The thickest shells from the largest mass-loss 
rates appear on the right-hand side of the diagram, where they exhibit a notably smaller 
latitude distribution. 
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